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Spatial interpolation of a surface electromyography (sEMG) signal from a 
set of signals recorded from a multi-electrode array is a challenge in 
biomedical signal processing. Consequently, it could be useful to increase 
the electrodes’ density in detecting the skeletal muscles' motor units under 
detection's vacancy. This paper used two types of spatial interpolation 
methods for estimation: Inverse distance weighted (IDW) and Kriging. 
Furthermore, a new technique is proposed using a modified nonlinearity 
formula based on IDW. A set of EMG signals recorded from the noninvasive 
multi-electrode grid from different types of subjects, sex, age, and type of 
muscles have been studied when muscles are under regular tension activity. 
A goodness of fit measure (R2) is used to evaluate the proposed technique. 
The interpolated signals are compared with the actual signals; the Goodness 
of fit measure's value is almost 99%, with a processing time of 100msec. 
The resulting technique is shown to be of high accuracy and matching of 


spatial interpolated signals to actual signals compared with IDW and 
Kriging techniques. 
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1. INTRODUCTION 

The surface electromyography (EMG) signal is the spatial and temporal summation of all active 
motor units’ action potential under regional detection’s vacancy [1], [2]. A high-spatial-precision EMG 
electrode, mounted on the skin surface of skeletal muscle, allows detection and estimation of the size of a 
single motor unit (MU). Moreover, analyzes neurological properties like conduct velocity single of the motor 
unit action potential (MUAP) activity when a muscular tissue is under contraction [3]. This paper’s main 
approach to enhance the surface EMG signal is based on spatial interpolation. The simple interpolation 
approach is to find the closest data value and to give the same value. This approach is not likely to be used as 
a linear interpolation [4]. Estimating missing data is an essential issue in signal processing [5]. Spatial 
interpolation has been commonly and frequently utilized in many analyses to obtain surface data dependent 
on a set of sampled points, like temperature, soil properties, and precipitation [6]. Consequently, spatial 
interpolation is one of the engineering tools to estimate spatial missing information in many signals like 
EMG. There are many researchers in both of those hot fields, enhancement of spatial detectors, and 
optimization of Spatial Interpolation, in the last few years. 

Abser et al. [7] illustrate a comparison between two interpolation methods, raw spectrum 
interpolation, and smooth spectrum interpolation, for error estimation signal-to-60Hz-noise ratio surface 
EMG signals. The results show that the smoothed spectrum is relatively better than the raw spectrum, but 
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both methods have estimation errors <10% relative to signal power. Hu et al. [8] proposes an interpolation 
method (RicFlow) of correspondence for noise cancellation. RicFlow has the smallest expected prediction 
error (EPE) relative to the other methods. The proposed method has (5.620) EPE. However, it still failed 
when low-density matches as input are applied. Martinez et al. [9] proposes a method to estimate grasp force 
(GF) using 16 EMG channels. Two Reduction methods, the Fix-EN and LassoG-EN, are used to select the 
number of features, to achieve between 90% and 99% of the initial number of features. According to the 
training windows (TW) length, the prediction improved until a stable performance was reached around 400 
ms after the grasp force (GF) onset. Li et al. [10] proposes an optimized method of inverse distance weighted 
spatial interpolation (IDW) called an adjusted inverse distance weighted spatial interpolation method 
(AIDW). The new method takes into account and position of sample points on the interpolated point due to a 
coefficient (K). The coefficient is added to set the sample point's distance weight due to its shielded effect in 
sample point positions. The proposed method (AIDW) has less mean related error (6.5%) than the IDW 
method (7.1%). Fan et al. [11] illustrate the evaluation effects of a uniform grid cell's size on the efficiency 
of k-nearest neighbor (kNN) search to locate neighbors fast. The standard deviation of points’ coordinates 
was used to measure the spatial distribution of scattered points. The relationships between the grid cell’s 
optimal size and the standard deviation of scattered points’ coordinates were fitted. This paper shows that 
with the increase of the standard deviation of points’ coordinates, the grid cell’s relatively optimal size 
decreases and converges. 


2. SPATIAL DETECTORS 

The surface EMG signal is the spatial and temporal summation of all active motor units’ action 
potential under regional detection vacancy. A high-spatial-precision EMG electrode mounted on the skin 
surface of skeletal muscle under contraction allows detecting and estimating the size of a single motor unit 
(MU), in addition to analyzes of neurological properties like conduct velocity of a single motor unit action 
potential. It is based on the multi-electrode; it consists of multiple electrodes arranged in a set of one or two 
dimensions, with a distant space between the electrodes of 4 millimeters [12]-[16]. Figure 1 shows a 
multi-electrode grid, designed, and developed in the electronic circuits laboratory\department of Electrical 
Engineering, University of Technology, Baghdad. It consists of (4x4) stainless steel pin-electrodes build on a 
plastic base and installed on a rubber layer to protect the skin of the subjects from damage. Because the 
signal’s amplitude may vary from less than 50 uV up to 30 mV, each picked signal is passed through a pre- 
amplifier, and the whole package is covered then by aluminum grounded shields to protect them from the 
noise. That was the detection, pre-amplification, and transduction stage of the signals. 


Figure 1. Electrode array with (4x4) electrodes 


An electrode-amplifier circuit that detects the sEMG signal amplifies and eliminates some of the 
unwanted noise. The electrode amplifier has a high input impedance to restrict current drawn from the 
subject and reduce signal leakage and attenuation. Without a change in output voltage, the low output 
impedance drives the following electronic circuit stage. The input impedance of modern operational 
amplifiers is over 100 MQ, and the output impedance is 100 Q. AD620 instrumentation amplifiers improve 
electronic efficiency (increase common-mode rejection ratio, decrease noise) and minimize inherent noise in 
subsequent stages. The amplifier is located very close to the electrodes, providing the most reliable sEMG 
information before amplification and reducing cable movement noise. Many researchers prefer this design 
because buffering the acquired signal near the source, amplifying the signal early in the process, and 
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minimizing the physical measurements of a system by only including a pre-amplification stage near the 
subject attenuates noise and motion artifacts. For avoiding the intrusion of hazardous current into the subject, 
all components’ power supply on the subject’s side of the circuitry is isolated from the ground. 

The second stage is to pass the 16-transduced analog signals to a 16-channel analog 
multiplexer/demultiplexer (CD74HC4067), multiplexer (CD74HC4067), which uses four digital pins to 
control the flow of one pin to 16 others. Next stage, the selected signal will be denoised by a band-pass filter 
(15-500 Hz) and then amplified. Finally, for the signal to be readable by the data acquisition system, a 
shifting level controlled by DC reference voltage (2.5 V) is applied. 

Data acquisition system that is used for the proposed approach is based on the Arduino Uno 
platform. Arduino Uno microcontroller module has a USB connection interface. The Arduino IDE and 
libraries are open sources, and there are a large variety of accessories, shields that drop onto the 
microcontroller. It has a variety of I/O pins, including analog, digital, and postmenopausal women (PWM) 
ports. The block diagram of the system is shown in Figure 2. 


Sensing and 16 channels of Band-pass 


Pre-Amplifying Analog switch Filter 


Arduino Shifting Level Amplifier 


Figure 2. The block diagram of the system 


When the skeletal muscle undergoes a voluntary contraction through nerve fiber, a possible time 
distribution variable of the potentials is detected on the skin's surface. Multi-electrode grids are used to 
observe the distribution of spatial potentials generated. From this, both localization of the potential source 
and its propagation in time can be obtained. The potentials detected by the grid's various electrodes can be 
considered the instantaneous spatial samples of all active MUs generated under the detection’s vacancy by 
activating the muscle fibers. This potential distribution contributes in various ways, such that the MUs 
located near the surface of the skin contribute to the potential output distribution higher potential [13]. 


3. SPATIAL INTERPOLATION 

Many interpolation techniques are available to predict the unknown area’s values like spline, 
Kriging, trend, natural neighbors, and inverse distance weighting. The biggest question among the 
researchers is to find out which is the desirable method for interpolating data. Several discussions and 
researches have been conducted to find out the best interpolation by Mahalingam et al. [17]. Here, we present 
the idea of a variable signal over time, which represents the instant spatial samples of the possible dispersion 
of the muscle fibers. Two widely and robust techniques are chosen to compare with our proposed technique. 


3.1. Inverse distance weighted 

Spatial interpolation techniques are used to determine the unknown places using known location 
values. This paper presents a modified nonlinear IDW interpolation technique to estimate a missing EMG 
signal from available neighbor’s signals. Normality and stationary assumptions are supposed to constrain the 
interpolation process within the data range. IDW is a robust and widely used technique to predict Olena’s 
unknown point [18]. It is based on the presumption that the interpolating area should be affected essentially 
by the proximal points and minimal of the far away points. The interpolating point is a weighted rate of the 
nearby points, the weight assigned to each nearby point decreases as the distance from the interpolation point 
to the neighbor points increases. It is calculated as (1) [18]. 


F(u) = >, wru a 
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Where, f(u) is the inverse distance estimation at location u, n is the number of locations, w;, 1=1, ...., n, are 
the weight attributed to each sample location, and f(u;) are the instantaneous data at the location i. The 
weights w; are determined by (2): 


1 
Ga) 
mp (t) 
i=1Gp) 
(2) 


where, (wi) are the euclidian distance between the interpolated data’s location and used sample points, p is the 
distance value's power. Observe that: 


n 
> w=1 


i=l (3) 


3.2. Kriging 

Like IDW, Kriging uses a weighting, which assigns more influence on the nearest data points in the 
interpolation of values for unknown locations. However, Kriging is not deterministic but extends the 
proximity weighting approach of IDW to include random components where the function does not know its 
exact location [19]. Kriging depends on spatial and statistical relationships to calculate the surface. The 
two-step process of Kriging begins with semi variance estimations [20]. i) variance=(zi-z;), ii) semi 
variance=variance/2, iii) the interpolation by using semivariogram. 

The advantages of this method are the incorporation of variable interdependence and the available 
error surface output. A disadvantage is that it requires substantially more computing and modeling time, and 
Kriging requires more input from the user (values and weights from every point) [21]. Kriging belongs to the 
family of linear least squares estimation algorithms. Kriging aims to estimate the value of an unknown 
real-valued function, E, at a point, (x;y), given the values of the function at some other points, 
fxl;y1);(42;y2);(43;y3)g. A Kriging estimator is said to be linear because the predicted value E(x;y) is a 
linear combination that for each pair of points separated by distance (h): i) take the difference between the 
attribute values, 11) square it, iii) add to sum, and iv) divide the result by the number of pairs. As written: 

y(h) = 


[Zi — Zj |? 


1 
N(h 
INO | Aton (4) 


4. SYSTEM SET-UP 

Four healthy subjects, males, and females (age 24-55 years) have different types of muscles, 
forearm flexor, extensor, Biceps Brachii, and m. Abductor Pollicisbrevis muscles are investigated. The study 
area consisted of a 3x3 dimensional arrangement of dry surface electrodes with a distant space of 5 mm. Each 
electrode (1-9) will measure the muscle activity (monopolar potentials), EMG1, to EMG9 under the area of 
sensing concerning ground, as shown in Figure 3. The detected nine signals are amplified and fed into a 
computer via an A/D converter. For each test, a period of 5-second signals is sampled at a 4 kHz rate and 
stored in the computer. A processing period of 100 msec is considered appropriate for EMG signal 
processing and applications, as mentioned by Harba and Ibrahim [22]. Figure 4 shows the recorded nine real 
signals for the subject two configurations of the spatial electrodes used. 


le +e 70 

T 5mm 
24 50 Be 
30 Se % 


Figure 3. The arrangement of the electrodes 
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Figure 4. Recorded nine real signals for subject 2 
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5. IDW AND KRIGING TECHNIQUES 

The study area of interest is to use the IDW and Kriging interpolation techniques to interpolate the 
instantaneous value of the signal (EMGS) in the area surrounding by electrode locations (1,7,3,9), from the 
instantaneous observation of the eight-neighbor signal’s (EMG1, EMG2, EMG3, EMG4, EMG6, EMG7, 
EMG8, and EMG9). A goodness of fit (R2) is applied to measure the discrepancy between the interpolated 
and the real values of the signal (EMG 5), 


0> R?>1 (5) 


where, 0: miss fitting, 1: good fitting. 

Using the IDW technique, the power (p) of the Euclidian distance wi shown in (2), was investigated. 
Eight different power values (0.5, 1, 2, 5, 10, 15, 25, and 30) are selected to find the suitable value of the 
instantaneous interpolated signal (EMG5). A standard exponential model Kriging method is chosen to 
interpolate the missing EMG signal at point 5 from points (1, 2, 3, 6,9,8,7, and 4) as shown in Figure 1, 
which needs to estimate the nugget, sill, and the range used to describe the semivariance for each set of 
signals and subjects. Table 1 shows the Goodness of fit for both methods, offering the best power of distance 
for IDW. Results show that a better fit For the IDW method could be achieved as the power increased and 
settled in the order of (p=15) for subject 1 (p=12) for subject 2 (p=8) for subject three, and (p=5) for 
subject 10. In comparison, the Kriging method shows a slight improvement than the other. Figure 5 shows 
the real and interpolated EMG signals with a total processing delay time of 100 msec., When the muscles are 
under moderate contraction, for IDW (p=5) and Kriging interpolation method. 


Table 1. IDW Interpolation for different module power 


Interpolation methods The Goodness of fit (R°) 
Sub 1 Sub 2 Sub 7 Sub 10 
IDW 0.8920 (p=15) 0.7835 (p=12) 0.8387 (p=8) 0.7808 (p=5) 
Kriging 0.9215 0.80137 0.8543 0.8001 
$ 
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Figure 5. The real and interpolated signals using IDW and Kriging methods for subject 7 


6. THE PROPOSED TECHNIQUES 

Spatial data interpolation works based on the first law of geography formulated by Waldo Tobler 
[23]. The law provides that “all the features on the earth’s surface have a relationship between one another, 
but the near features are having more relation than features away.” However, for sEMG, that assumption is 
either valid if the detection area belongs to a single motor unit’s activity or false if the detection area belongs 
to more than one motor unit. This assumption depends on the motor unit's size, which differs from muscle to 
muscle, as mentioned by Nandedkar et al. [24], and the detector’s surface area. For this reason, further 
measures have been taken to study the effect of the presence of more than one motor unit under the area of 
detection. After evaluating the previous results, the number of recorded signals surrounding point 5 may have 
a weak effect on interpolating the missed signal, as shown in Figure 4. For that reason, two opposite spatial 
interpolation locations are implemented using IDW; Table 2 shows the effect of two opposite points to 
interpolate the signal at the location (5). The Goodness of fit is used to evaluate the interpolated signals 
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concerning the real one. Results show that a better fit is achieved for locations 4 and 6 for all subjects, and 
that is true when the electrode axis is mounted parallel to the muscle’s fibers, but this may not be valid if it is 
out of axis. In other cases, the fitness gets worse when the electrode axis is diverted toward the muscle fibers. 
There may be more than one motor unit under the area of detection. These results are seen to be correct as 
Figure 6 the amplitude of real and interpolated with order (p=8) for the sEMG signals’ cases, Figure 6(a) 
locations 1 and 9, Figure 6(b) locations 2 and 8, Figure 6(c) locations 3 and 7 and finally Figure 6(d) for 
locations 4 and 6, with a total processing period of 100 msec, under a moderate muscle contraction. The best 
fit is obtained for interpolating the instantaneous signal at location 5 from instantaneous signals at locations 4 
and 6, and the worst fit is obtained at locations 3 and 7, while the other gives a moderate fit. That explains 
that there is more than one source of activity under the vacancy of detection, which will affect the 
instantaneous value of the interpolated signal. 

From (1), the distance is only the criteria used to interpolate the signal, which satisfies the basis of 
the first law of geography formulated by Ouabo et al. [20]. In contrast, our results show that the detection 
area may be affected by more than one activity source for online interpolated signals. It was taken into 
consideration to develop a robust approach to interpolating the missed signal. For this basis, distance, and 
Goodness of fit (R2) between the real and interpolated signals were taken into consideration; a nonlinearity 
approach to the inverse distance is applied. A coefficient (G), this nonlinearity is present, and the (1) is 
modified to: 


oe Kw Gif (ur) (6) 


Gi has a coefficient derived from the normalized to one; the Goodness of fit is shown in Table 3. 

A gain (K) is made to overcome the interpolated signal's weakness by multiplying by the 
coefficients (Gi), normalized to one, to obtain the highest efficiency possible. The gain factor (K) is obtained 
by using the simple linear regression eguation stated by [25], [26]: 


Y=a+bx (7) 


The y“intercept is represented by (a), and the slope of the line is indicated by b. We obtain the simple linear 
regression equation using site sink source profiles (SSSP); by entering the x values as the quantized absolute 
mean value of each subject’s eight surrounding signals. The output simple linear regression is (8): 


9 = 4.985 + 2.685 x (8) 


Which is more suitable for the best signal interpolated gain (K) for all cases. Table 4 illustrates the Goodness 
of fit between the interpolated signals and real signal for the IDW and the modified IDW (proposed 
coefficient G). The result shows that the modified IDW spatial interpolation method has better-predicted 
accuracy than the standard IDW spatial interpolation, with an enhancement of (9.85-26%). 

The time responses, shown in Figure 7, show the power of matching between the spatial 
interpolation and the real sEMG signal, despite small distinction and the shape of the interpolated signal. The 
proposed technique shows reliability, confidentiality and can be applied to other types of biomedical signals. 
The non-normality shape of the signal distribution function and belongs to the non-convergence property is 
derived from the results of the conventional method in the case of non-contaminated distributions. 


Table 2. IDW for different points 


Subj. Points 4-6 3-7 2-8 1-9 
1 R? 0.9671 0.7027 0.7603 0.7945 
2 R? 0.9929 0.0650 0.7463 0.6998 
7 R? 0.9756 0.6405 0.6019 0.6584 
10 R? 0.9811 0.6450 0.7619 0.6230 


Table 3. The normalized Goodness of fit for Table 4. The Goodness of fit between the real interpolated 


different points (G) signals, before and after using the nonlinear factor (G) 
Subj. 4-6 3-7 2-8 1-9 Subj. _R° (before) __R? (after) Improvement  K 
1 0.9671 0.7027 0.7603 0.7945 1 0.8920 0.9803 9.8% 5.5 
2 0.9929 0.0650 0.7463 0.6998 2) 0.7835 0.9900 26.6% 5.4 
7 0.9756 0.6405 0.6019 0.6584 7 0.8387 0.9620 14.7% Su 
10 0.9811 0.6450 0.7619 0.6230 10 0.7808 0.9861 26.3% 5.9 
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Figure 6. The real and interpolated signals for cases, (a) 1-9, (b) 2-8, (c) 3-7, and (d) 4-6 
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Figure 7. The real and spatial interpolated signals for subject 7 


7. CONCLUSION 

A new nonlinearities spatial interpolation formula is applied to interpolate the instantaneous value of 
the sEMG signal located at a specified location, by using the instantaneous observation of the surrounded 
neighbor’s signals, to increase the electrodes’ density in detecting the skeletal muscles’ motor units under 
detection’s vacancy. It must be mentioned that the results obtained are contrary to our hopes that the resultant 
technique is achieved better performance than the simple IDW. The time responses, shown in Figure 7, show 
the power of matching between the spatial interpolation and the real sEMG signal, despite the small 
distinction and the interpolated signal’s shape. It can be concluded that the proposed technique shows 
reliability, confidentiality and can be applied to other types of biomedical signals. 
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Hardware limitation of electrodes like technical problems and cost is an important issue that faces 
the EMG surface mounted electrode industry. The limited aperture between electrodes illustrates the 
technical problems. In contrast, the proposed system is a good solution for both issues (electrodes limitation 
and cost), by our modified technigue to interpolate the desired signal extracted from known signals under the 
vacancy of the detection area. In other words, increase the precious detection of the electrodes. In the case of 
non-contaminated distributions, the suggested formula, which deals with the non-normality shape of the 
signal distribution function and belongs to the non-convergence property, derives from the results of the 
traditional method, so that the results are better than those of the conventional IDW method when dealing 
with contaminated distributions. The suggested formula has a high degree of accuracy between real and 
spatial interpolated signals regardless of different subjects or muscles, with an improvement in the Goodness 
of fit (9.85-26%). One of the significant recommendations is to reconsider the characteristics of the missing 
signals, of the studied cases, through surrounded neighbor signals’ instantaneous values. Furthermore, it is 
recommended that the suggested formula can be applied to the physical body’s characteristics, such as the 
study of the motor unit action potential. 
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